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Abstract 

It is known that the parity of reflection amplitude can either be even or odd under the mirror 
operation. Up to now, all the parities of reflection amplitude in the one-mode energy region are even 
under the mirror operation. In this paper, we give an example of odd parity for Andreev reflection 
(AR) in a three-terminal graphene-supercondutor hybrid systems. We found that the parity is even 
for the Andreev retroreflection (ARR) and odd for specular Andreev reflection (SAR). We attribute 
this remarkable phenomenon to the distinct topology of the band structure of graphene and the 
specular Andreev reflection involving two energy bands with different parity symmetry. As a result 
of odd parity of SAR, the SAR probability of a four-terminal system with two superconducting 
leads (two reflection interfaces) can be zero even when the system is asymmetric due to the quantum 
interference of two ARs. 
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Since the experimental realization of graphene jji], it has become an exciting arena for 
theoretical and technological investigations. [4] A number of new phenomena have been pre- 
dicted and verified experimentally. For instance, in the presence of magnetic field, it exhibits 
a distinctive half-integer quantum Ha., effect^ Its quasi-particles obey the Duae-hhe equa- 
tion and have relativistic-like behaviors. |4J Due to the relativistic effect, the Klein tunneling 
occurs where an incident electron in graphene can pass through a potential barrier with 
probability one.[5| Then a graphenep-n junction can be used to focus Dirac electron current 
with a negative refractive index. Q] 

Since good contacts between superconducting leads and graphene have been re- 
alized experimentally, 8j] the transport study through graphene based normal-metal- 
superconductor (GNS) heteroj unction becomes feasible. In the presence of a normal metal 
(graphene)-superconducting interface, an incoming electron converts into a hole and a cooper 
pair is formed that enters the superconductor. Due to the relativistic nature of the electron 
in graphene, the electron-hole conversion can either be intraband (within conduction or va- 
lence band) or interband (between conduction and valence bands). When the electron-hole 
conversion is intraband, it corresponds to the usual Andreev reflection (AR) 9] or Andreev 
retroreflection (ARR) because the reflected hole is along the incident direction. This ARR 
occurs for both relativistic and non-relativistic electrons. When the electron-hole conversion 
is interband, the reflected hole is along specular direction and a specular Andreev reflection 
(SAR) takes place, 1(| which can lead to novel phenomena as we will discuss below. 

It is known that the parity is a fundamental quantity in physics and reflection is a 
general physical phenomenon in nature. In this paper, we discuss the parity of reflection 
amplitude for graphene in contact with superconductor leads. In general, the parity of a 
reflection amplitude can be either even or odd when the system is under mirror operation. 
However, for all previous known reflection events, the reflection amplitudes in the one-mode 
energy region have even parity under the mirror operation. It is yet to find an odd-parity 
reflection event. In this paper, we found, for the first time, that the SAR amplitude has 
an odd parity under the mirror operation for zigzag graphene ribbons with even number 
of chains. This means that the phases of SAR amplitude for a graphene-superconductor 
hybrid system and its mirror system differ by ir. We attribute this phenomenon to the 
unique band structure of the graphene. Obviously this phase difference does not affect 
any observable quantities for each system. When two systems couple together, however, 
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this 7r phase manifest through quantum interference between two SARs. So this ir phase 
shift has important consequences for a four terminal device with two superconducting leads 
(see Fig.2(a)). When two superconducting leads are symmetrically attached to the device, 
the quantum interference of the left and right SAR leads to a destructive or constructive 
interference depending on whether the phase difference of superconducting leads is zero or 
7i. Importantly, when two superconducting leads are asymmetrically attached to the device, 
the same interference pattern occurs provided that the Dirac point Eq is in line with the 
condensate of superconducting lead. The quantum interference between pairs of the AR can 
be tuned by shifting the Dirac point, the asymmetry of the two superconducting leads, as 
well as the phase between two superconducting leads. Due to the odd parity of SAR, the 
interference pattern for SAR is phase contrasted to that of ARR where the parity is even. 

Before doing numerical calculation, we first prove that the phases of SAR amplitude of 
two systems (i) and (ii) in Fig. 1(a) differ by tt, i.e., the parity of SAR is odd under mirror 
operation. Note that for graphene systems electrons in valence and conduction band are 
usually referred as electrons and holes, respectively. In the presence of superconducting 
lead the reference point of electrons and holes is the Fermi level in the superconducting 
lead. In the following, we will refer electrons (holes) as electrons above (below) Fermi 
level in superconducting lead. Denote ip£ {ipt) the wavefunction of electrons in conduction 
(valence) band moving in +y direction and ip~ in -y direction in the zigzag graphene 
nanoribbon lead. It was known that under reflection P : x — > —x, ipf is symmetric while 
ipf- is anti-symmetric if the energy of electron is in the first transmission channel [ill] [see 
Fig.l(a)], i.e., 

P^{x, y ) = -^(-x,y). (1) 

which is one of the unique features of zigzag edge nanoribbons with even number of chains. 
Assuming the incident electron from the terminal-1, the wavefunctions for SAR in zigzag 
nanoribbon lead 1 or 3 of the system (i) can be written as 

4 ] = W e + + (2) 
where r n is the normal reflection amplitude, ti 3 is the transmission amplitude, r n A and 



ri3A are the Andreev reflection amplitudes with the reflected hole to the terminal- 1 and 3, 
respectively. Similarly the wavefunctions for the system (ii) are given by 

4 U) = We + +n 3 A^ (3) 

Since the system (i) is related to (ii) by the reflection operator P, we have ipa — Pipa with 
a = 1,3. Note that for SAR, the electron is in the conduction band while the hole is in the 
valence band, i.e., ip e = ip c and iph = ip v . From this relation together with Eqs.(DQ), 02]), and 
02]), we obtain 

ruA = -fiiA, nsA = -ri3A 

ru = f n , t 13 = ti 3 (4) 

Note that the origin of this 7r phase shift (odd parity) is the interband conversion from the 
electron to the hole. Therefore the ir phase shift does not occur for ARR since it involves only 
intraband conversion. Now we verify this statement numerically using a tight-binding model 
(see below for detailed description of the model and numerical procedure). The numerical 
results of AR probability -Riiyi(i3yi) = |?"iiA(i3A) | 2 for two systems are shown in Fig.l(b). As 
expected the AR probability are exactly the same for two systems. However, the phase 
of AR amplitudes rnA(i3A) denoted as 3>iiq 3 ) are different. It is shown in Fig.l(c) and 
Fig. 1(d) that ARR amplitudes {\Eq\ > \Ep\, with \Ep\ = 0.5) are the same for two systems 
in Fig.l(a) while the SAR amplitudes {\Eq\ < \Ep\) have a ir phase shift. It confirms the 
odd parity for interband electron-hole conversion, which comes from the distinct topology 
of the band structure of graphene. 

To see the consequence of the odd parity of SAR, we examine a symmetric four-terminal 
device with two superconducting leads depicted in Fig.2(a) (by setting asymmetry SN = 
and phase difference 5(f) = 0). For this system, two beams from terminal-1 has a ir phase 
shift due to odd parity of SAR and interferes destructively at terminal-3 giving rise to a 
vanishing SAR coefficient. However, we can arrive the same conclusion using symmetry 
argument as follows. Since the system is symmetric with respect to x = 0, we must have 
^13,4 = ^i3A when the reflection operation along x-direction is applied. While from Eq.(4), 
r 13 A = — ri3A- So the AR probability R13A = |^i3^| 2 for SAR can also be zero from symmetry 



point of view. 



121 ] Therefore we conclude that the symmetric device can not be used to test 



the odd parity of SAR. In the following, we demonstrate that due to the ir phase shift the 
destructive interference still occurs in a four-probe devices with two superconducting leads 
attached asymmetrically and hence can be used to test the odd parity of SAR. 

For this purpose, we consider an asymmetric four-terminal device consisting of a zigzag 
graphene ribbon with two superconducting leads as shown in Fig.2(a). The Hamiltonian of 
the graphene is|l3!] H = tia\ a\ — X]<ij>^ a i a j- Here a\ and a\ are the annihilation and 
creation operators at site i, ei is the on-site energy which can be controlled experimentally 
by the gate voltage^, and the hopping constant t = 2.75eV represents the nearest carbon 
bond energy The pair potential (energy gap) of superconducting terminal-/^ with (3 = 2,4 
is = A^e^ with A 2 = A 4 = A ~ lmeV. In numerical calculations, 12| we fix Fermi 
energy Ep and tune the Dirac point E . We have used A as the energy unit. 

Now we study the interference between two ARs from GNS junctions as shown in Fig. 2(a) 
in which two superconducting leads 2 and 4 are asymmetrically attached to the zigzag 
nanoribbon. The horizontal distance SN between two GNS junctions measures the asym- 
metry of two GNS junctions. The scattering process can be qualitatively understood as 
follows. For simplicity, we assume 02 = 04 for the moment. As shown schematically in 
Fig.2(a), for SAR the particle-like electrons in terminal-1 split into two beams and are scat- 
tered separately by two GNS junctions (green horizontal lines) as holes that finally recombine 
at terminal-3. We examine the total phase accumulated for each beam that involves the fol- 
lowing three processes. Before reaching the first GNS junction (denoted by the left vertical 
green line) two beams of electrons propagate with the same momentum k x . After reaching 
the second GNS junction (denoted by the right vertical green line) two beams of holes also 
propagate with the same momentum k' x . Obviously phases accumulated in the above two 
processes for both beams are the same. Between them two beams propagate with different 
momenta k x and k' x . Hence the phase difference between two beams is = (k x — k x )Sx with 
Sx = bSN, where b = a/3o and a the lattice constant. This phase difference can be tuned 
by varying the Dirac point E or the asymmetry SN giving rise to a complicated interfer- 
ence pattern (see Fig.2). In particular, this phase difference can be zero if (k x — k' x ) = 
(i.e.,Eo = 0) or SN = 0. In general, the total phase difference is = (k x — k x )5x + 02 — 04- 

Interference pattern of AR probability R\za for system depicted in Fig. 2 (a) with pair 
potential phase difference of two superconductors Sep = and tt (Sip = (f2 — tp^) are then 
plotted in Fig.2 (b) and (c), respectively. For Fig.2 (b) following observations are in order: 
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(1) For the geometrically symmetric system (SN = 0), the interference is always destructive 
with zero Riza as long as |_Eb| < |-£ i f|.[12| Clearly this is due to the n phase shift depicted 
in Fig. 1(d) and is consistent with the band selection rule, 11] (2) When Dirac point Eq is 
in line with the condensate energy of the superconductor, i.e., when E = 0, -R13A is again 
zero no matter what value SN assumes. This means that there is a completely destructive 
interference between two beams scattered by two GNS junctions attached asymmetrically to 
the graphene nano-ribbon. This behavior can be understood as follows. When Eq = the 
incoming electron and reflected hole have the same propagating momentum k x and thus path 
1 and 2 in Fig. 2 (a) experience the same quantum phase k x Sx except at the superconducting 
leads. Hence the total phase difference is only due to the 7r phase shift between two SARs. 

(3) R13A is an even function of Dirac point Eq because of the electron-hole symmetry in 
graphene. Due to the geometric symmetry, R13A is also an even function of asymmetry SN. 

(4) For nonzero Ep, the closer the Dirac point Eq to Ep, the more rapidly -R13A oscillates 
as we vary SN. This is because the difference of propagating momentum k x — k' x increases 
monotonically as Eq approaches to Ep. (5) When Eq is in the vicinity of Ep, R13A can 
reach 0.9 which is much larger than that when \Eq\ > \Ep\. This is because when Ep is 
very close to Eq, the edge states of zigzag ribbon begin to contribute, then electron is easier 
to be scattered by two GNS junctions located also at edges of zigzag ribbon. Considering 
the pseudo-spin conservation, large R13A is always found in the region of \E \ < \E F \, i.e., 
the SAR region. (6) There is an overall fine oscillation with a period of SN = 3b. Similar 
behavior was also found in zigzag ribbons with a p-n junction where the conductance is 
determined by the relative displacement S along the p-n junction. jl4] In Fig.2(c) with the 
superconducting phase difference Sip = tt, we see that the interference pattern is contrary to 
Stp = [Fig.2(b)] where the constructive interference becomes destructive and vice versa. 

To further analyze the interference pattern, we plot in Fig. 3 (a) the total R13A vs Dirac 
point Eq for different asymmetry SN with the phase difference between two superconducting 
leads Sip = [main panel of Fig.3(a)] or S<p = tt [inset of Fig. 3(a)]. Clearly the interference 
(oscillatory) pattern occurs only for asymmetric systems (SN 7^ 0) with oscillation frequency 
proportional to SN. When pair potential phase difference Sip = tt is introduced, the inter- 
ference pattern reverses, and -R13A with SN = becomes the envelop function of R13A for 
all nonzero SN. In Fig.3(b) we plot -R13A vs SN for different widths W of nanoribbon. It 
is shown clearly that -R13A is a periodic function of SN with larger periodicity for larger W. 
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In the inset of Fig.3(b) we plot this period versus the width for different Eq. The period 
P is obtained in two ways: (1). from the expression P = 2ii/(k x — k' x ) where the momenta 
k x and k' x can be obtained from the band structure for a given E (black symbols). (2). 
directly from main panel of Fig.3(b) (red solid circle). From the inset, it clearly shows that 
two periods are exactly the same giving strong evidence that the interference pattern of AR 
probability are indeed from two reflected hole beams. 

Finally, the interference pattern of AR probability Rua is also studied (not shown). We 
found that only ARR probability -Rua {\Eq\ > Ef = 0.2A) exhibits interference pattern. We 
note that since there is no tc phase shift involved in ARR, when 5N = reflected electrons 
through two GNS junctions interfere constructively when 5tp = and destructively when 
Sip = 7T which is in contrast to SAR in Fig. 2. In fact, interference patterns of SAR and ARR 
are always phase contrast not only for SN = but also for all other SN. 

To test the odd parity of SAR experimentally, it relies on the fabrication of high quality 
zigzag graphene nanoribbons. It has been achieved by several laboratories using differ- 
ent methods last year including the method to unzip the multi-walled carbon nanotube 
(CNT),[l^] the anisotropic etching by thermally activated nickel nanoparticles, 16] and use 



reconstruction of the edge to make zigzag graphene nanoribbons. [17J In view of the above 
experimental breakthrough, we expect that the setup to test our predicted phenomenon can 
be realized experimentally. 

To reduce the experimental challenge, we have considered an unzipped CNT device, i.e., 
(n,n) CNT-zigzag graphene- (n,n) CNT, obtained by unzipping a few unit cells in the central 
part of an armchair CNT which has been achieved experimentally. 15f] For this system, the 
wavefunction in the armchair CNT has the same symmetry as that of the zigzag graphene 
ribbon. Following the same procedure leading to Eq.(4), we have shown that the unzipped 
CNT in contact with a superconducting lead has the odd parity under mirror operation. 
Similar conclusions drawn from GNS can be obtained for unzipped CNT with two super- 
conducting leads. 

In conclusion, up to now, the parity of reflection amplitude was found to be even under 
the mirror operation. Here we have provided an example of odd parity for the reflection 
amplitude, the SAR amplitude in the zigzag graphene-superconductor hybrid system. This 
odd parity is due to the combination of unique band structure of the graphene and the 
electron-hole conversion involving two energy bands with different parity symmetry. The 
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signature of odd parity of SAR can be found from the quantum constructive interference in a 
four terminal system with two superconducting leads attached asymmetrically. Furthermore, 
the interference pattern due to odd parity of SAR is phase contrasted to that of ARR where 
the parity is even. 

Acknowledgement We gratefully acknowledge the financial support from a RGC 
grant (HKU 705409P) from the Government of HKSAR and from NSF-China under Grant 
Nos.10974236 and 10821403. 



[1] *e-mail: jianwang@hkusua.hku.hk. 
[2] ^e-mail: sunqf@aphy.iphy.ac.cn. 

[3] K. S. Novoselov et al, Science 306, 666 (2004); K. S. Novoselov et al, Nature (London) 438, 

197 (2005); Y. Zhang et al, Nature (London) 438, 201 (2005). 
[4] C. W. J. Beenakker, Rev. Mod. Phys. 80, 1337 (2008); A.H. Castro Neto et al, Rev. Mod. 

Phys. 81, 109 (2009); A. Rycerz et al, Nature Phys. 3, 172 (2007). 
[5] M. I. Katsnelson et al, Nature Phys. 2, 620 (2006); A.-F. Young et al, Nature Phys. 5, 222 

(2009); A. V. Shytov et al, Phys. Rev. Lett. 101, 156804 (2008). 
[6] V. V. Cheianov et al, Science 315, 1252 (2007). 
[7] Y. Xing, J. Wang and Q.-F. Sun, Phys. Rev. B 81, 165425 (2010). 
[8] H. B. Heersche et al, Nature 446, 56 (2007); F. Miao et al, Science 317, 1530 (2007). 
[9] A. F. Andreev, Sov. Phys. JETP 19, 1228 (1964). 
[10] C. W. J. Beenakker, Phys. Rev. Lett, bf 97, 067007 (2006). 
[11] J. Nakabayashi et al, Phys. Rev. Lett 102, 066803 (2009). 

[12] S. Cheng et al, Phys. Rev. Lett. 103, 167003 (2009); Q.-F. Sun and X.C. Xie, J. Phys.: 

Condens. Matter 21, 344204 (2009). 
[13] D. N. Sheng et al, Phys. Rev. B 73, 233406 (2006); Z. Qiao and J. Wang, Nanotechnology 

18, 435402 (2007). W. Long et al, Phys. Rev. Lett. 101, 166806 (2008); J. Li and S.-Q. Shen, 

Phys. Rev. B 78, 205308 (2008). 
[14] A. R. Akhmerov et al, Phys. Rev. B 77, 205416 (2008). 
[15] L. Jiao et al, Nature (London), 458, 877 (2009). 



8 




FIG. 1: (Color online) Panel (a): zigzag ribbons with even number of chains (gray honeycomb) 
attached by a superconducting lead on the left and right (orange honey comb), respectively. For 
SAR the incoming electrons (red arrow) are scattered by the GNS junction (green solid line) as holes 
(blue arrow). The corresponding wave functions at sublattice "A" (solid circle) and "B" (hollow 
circle) for the lowest subband in conduction band (bottom) and the highest subband in valence 
band (top) are shown schematically. Panel (b): AR probability from terminal-1 to terminal-1 Rua 
and to terminal-3 R\za vs. Dirac point Eq. Panels (c) and (d): AR phase (c) and (d) of 
two systems in panel (a) and their phase different — ^11(13) vs - Eo- 

[16] L. C. Campos et al, Nano. Lett. 9, 2600 (2009). 
[17] C- O. Girit et al, ibid, 323 1705 (2009). 
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FIG. 2: (Color online) Panel (a): sketch of AR interferometer in which the zigzag ribbon is 
asymmetrically attached by two superconductor lead-2 and 4. Electrons in terminal-1 can be 
Andreev reflected into terminal-3 by either top or bottom GNS junction (horizontal green lines). 
Panel(b) and (c): the contour plot of R13A vs Dirac point £0 and asymmetry 5N. The phase 
difference of two superconductor leads 8<p is zero in panel (b) and tt in panel (c). The other 
parameters: Fermi energy Ef = 0.8, number of chains in zigzag ribbon N = 40 corresponding to 
width 60a, the width of superconductor lead Ws = 106, where b = y/?>a. 



10 




-0.6 -0.3 0.0 0.3 0.6 




asymmetry 5N (b) 

FIG. 3: (Color online) Panel (a): With fixed Fermi level Ep = 0.8, total AR probability R\%a vs 
Dirac point Eq for different asymmetry SN. In the main panel Sip = 0, and while Sip = tt in the 
inset. Panel (b): R\?,a vs asymmetry SN with Eq = 0.3i for different width W from 10 x 3a to 
38 x 3a with the interval 2 x 3a along the black arrow. Inset panel: the hollow signs are the period 
P obtained from the main panel and the solid red circles are the period P from the energy band 
with the expression P = 2n/(k x — k' x ). The other parameters: Sep = 0, Ep = 0.8. 
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